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Unconventional spin fluctuations in the hexagonal antiferromagnet YMnOa 
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We used inelastic neutron scattering to show that well below its Neel temperature, Tn, the two- 
dimensional (2D) XY nearly-triangular antiferromagnet YMnOs has a prominent central peak asso- 
ciated with 2D antiferromagnetic fluctuations with a characteristic life time of 0.55(5) ps, coexisting 
with the conventional long-lived spin-waves. Existence of the two time scales suggests competition 
between the Neel phase favored by weak interplane interactions, and the Kosterlitz-Thouless phase 
intrinsic to the 2D XY spin system. 

PACS numbers: 75.30.Ds, 75.40.Gb, 75.25.+Z, 75.50.Ee 
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Geometrical frustration and low dimensionality are the 
two key concepts in the statistical physics that provide 
unusual spin dynamics as well as phase transitionsi*^. 
The simplest realization of the two concepts is two- 
dimensional triangular lattice antiferromagnets (2DT- 
LAFM). A particular interest is placed on the XY spin 
system (2DXFTLAFM), where its ground-state mani- 
fold has the continuous degeneracy associated with U{1) 
global spin rotations, as well as the discrete Ising-like de- 
generacy due to Z2 chirality configurations^. Because of 
the U{1) symmetry, the well-known Kosterlitz-Thouless 
(KT) phase involving vortex binding^ is expected at low 
temperatures. Experimentally, little evidence can be 
found in the literature due to the lack of good model 
systems^.. 

Hexagonal rare-earth manganites iJMnOa {R =Y, Lu, 
and Sc; space group PQ^cm) can be good candidates for 
the 2DTLAFMS. The Mn^+ {S = 2) ions form nearly 
triangular networks in 2: = and 1/2 layers, stacking 
with the ABAB sequence (see Fig. The layers 

are well separated by R and O ions, suggesting good 
2D character in the ab plane. Their bulk susceptibil- 
ity data show that despite strong antiferromagnetic in- 
teractions the magnetic ions order long-range at much 
lower temperatures Tn than the magnetic energy scale 
inferred by the Curie- Weiss temperatures 0cw (for in- 
stance, Qcw = -705 K and Tn = 70 K for YMnOa)^^*^. 
Previous powder neutron diffraction studies showed that 
the spins at the lowest temperature formed the so-called 
120° structure in the a6-plane coinciding with the ground 
state for 2DTLAFMs, and the frozen moments, (M), 
were reduced from the expected value for the fully polar- 
ized Mn3+ {e.g., (M) = 2.90(2)^b < gS^s for YMnOg 
{S = 2))2iiia. The reduction in Tn and (M) is a 
signature of strong spin fluctuations due to geometri- 
cal frustration and/or low dimensionality in the sys- 
tems. A broad peak was additionally observed at finite 
Q around Tn, indicative of strong short-range spin cor- 
relationsSiifl. However, due to intrinsic limitations of the 
powder diffraction technique, further experimental stud- 
ies, especially inelastic single crystal neutron scattering 




FIG. 1: Schematic drawing of Mn positions and spin ordering 
in YMnOs. Filled (open) circles represent Mn positions in the 
z = {z = 1/2) plane, whereas dotted parallelogram shows 
the magnetic unit cell that is identical to the chemical unit 
cell. Lattice constants are a = 6.140 A and c = 11.393 A. The 
lattice is weakly trimerized with the intra-trimer and inter- 
trimer Mn-Mn distances of 3.42 A and 3.62 A, respectively^, 
which is exaggerated in drawing. 



measurements are necessary to understand the nature of 
the spin excitations. 

In this paper, we report for the first time inelastic 
neutron scattering measurements on powder and single- 
crystal samples of the hexagonal rare-earth manganite 
YMnOa. We have found that YMnOa is a good model 
system for the 2DAyTLAFM with weak trimerization. 
Our most important finding is that in the Neel phase 
there are fast 2D spin fluctuations with a characteristic 
time scale of 0.55(5) ps in addition to the conventional 
long-lived spin wave excitations. The inverse of dynamic 
correlation length associated with the fast 2D spin fluc- 
tuations has similar T-dependence as that expected for 
the KT phase in a 2D XY spin system, suggesting that 
the spin fluctuations are reminiscence of the KT phase. 
The coexistence of the long-lived magnons and the fast 
2D spin fluctuations also suggests competition between 
the Neel phase and the KT phase in this quasi-2D XY 
spin system. 

A 50 g powder sample and a 2 g (05 mm x 22 mm) 
single crystal of YMnOs were used in our neutron scat- 
tering measurements. Methods of sample preparation 
were reported elsewhere^i. Neutron scattering measure- 
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FIG. 2: Color contour maps of the powder-averaged neutron 
scattering intensity versus magnitude of wave vector trans- 
fer Q and energy transfer hui at tiiree different temperatures 
spanning the phase transition at Tn = 70 K. 



ments were performed at the NIST Center for Neutron 
Research. Powder experiments were performed at the 
Disk Chopper time-of-flight Spectrometer (DCS) using 
an incident energy of Ei ~ 15.46 meV and single crystal 
experiments at the cold neutron triple-axis spectrometer 
SPINS and the thermal neutron triple-axis spectrome- 
ter BT9. At SPINS, pyrolytic graphite (PG) 002 reflec- 
tions were used for monochromator and analyzer, and 
a cooled Be filter was placed after the sample to elim- 
inate higher-order contamination. We used horizontal 
collimations of 80'-80' and a final energy i?f = 5 meV 
for most scans, while = 2.6 meV and 80'-40' were 
used when better energy resolution was needed. At BT9, 
the PG monochromator and analyzer were used with 
Ef = 14.7 meV, and a PG filter was used to get rid 
of higher order contamination. 

Fig. 2 provides an overview of the inelastic neutron 
scattering intensity I(Q,lu) for the powder sample at 
three temperatures. The powder averaged scattering in- 
tensity is related to the dynamic spin correlation function 
S''^{Q,uj) asii. 
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where F{Q) is the magnetic form factor for Mn'^+. For 
T > Tn, there is a cooperative paramagnetic continuum 
centered at Q = 1.2 due to fluctuations of small 
AFM clusters, as is commonly found in geometrically 
frustrated AFMe^^. By integrating I{Q,uj) over fiui and 
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FIG. 3: (a-d) Constant Q = (1,0,0) scans at four differ- 
ent temperatures, (e) Constant-Q scans at Q = (1,0,2) and 
Q = (1,0,1) at T = 7 K. The scan at Q = (1,0,1) was 
taken under the higher-energy-resolution configuration with 
Et = 2.6 meV, and its unit is arbitrary. Solid lines are fits to 
Eq. Q, whereas dashed lines represent the quasielastic part 
(see the text). 



Q, we obtained the sum rule of S{S + 1) = 5.2(5) at 
80 K, which is close to the expected value for dynamic 
Mn'^+(S' = 2) ions. This and the Q dependencei^ tell us 
that the scattering is magnetic. For T < Tn, as the mag- 
netic long range order develops, spectral weight at low 
energies gradually shifts to higher energies. At T = 4 K 
there is strong scattering above tiuj ~ 5 meV and weak 
scattering below. 

Next, to obtain Q-directional dependence of the mag- 
netic excitations, we have performed single crystal in- 
elastic scattering experiments. Fig. 3 shows the repre- 
sentative constant-Q scans at the antiferromagnetic zone 
center F, namely, Q — (1,0,0) and equivalent positions. 
For r > Tn, the cooperative paramagnetic continuum 
appears as a quasielastic peak centered at hu) = meV. 
For T < Tn, the quasielastic peak intensity decreases and 
two prominent magnon peaks develop at nonzero ener- 
gies. The energy values of the magnon peaks increase as 
T decreases, becoming hu = 2.3 and 5.3 meV at 7 K. 
A constant Q = (1,0, 1) scan with a better energy res- 
olution revealed an additional mode at Tiui — 0.22 mcV 
(Fig. 3(e)). 

We analyzed the observed spectra using the following 
scattering function with Lorentzians for the quasielastic 
and magnon peaks: 
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FIG. 4: Spin wave dispersion relations along {h, 0, 0), {h, h, 0) 
and (0,0,1). Circles and triangles represent the experimen- 
tally determined peak positions. (Different marks distinguish 
branches appearing from different magnetic Bragg positions.) 
Solid lines are the model dispersion relations explained in the 
text. Inset: schematic drawing of scan directions in the 2D 
plane. 



where — [1— exp(— Tiw/feBT)]"^. This function 

was convoluted with the instrumental resolution function 
to fit the observed spectra. 

Let us first discuss the magnon dispersion relations at 
T = 7 K << Tn. Fig. 4 shows the dispersion relations 
along a few high symmetry directions, obtained from sev- 
eral constant-Q and constant-?iw scans. To explain the 
observed dispersion relations, we introduce the following 
model spin Hamiltonian: 

^ = - E J^JS^■S,~DlY,iS^f-D2J2iS^■n^)^ (3) 

<ij> i i 

which consists of two inplane (Ji and J2) and two inter- 
plane (J{ and J2) interactions, and the easy-plane (-Di) 
and inplane easy-axis (£'2) anisotropics (see Fig. 1 for the 
definition of the interactions). The anisotropy D2, paral- 
lel to the spin directions (n^ — {Si) /\{Si)\), is necessary 
to reproduce the small (0.22 meV) gap at the antifer- 
romagnetic zone center, and is presumably due to the 
local structural distortion around Mn'^+. The conditions 
J( > and J'l > J2 are necessary for the particular in- 
terplane stacking in YMnOs to be the ground state. 

The model Hamiltonian is linearized using the 
Holstein-Primakoff approximation, and numerically diag- 
onalized to obtain one-magnon dispersion relations using 
the standard equation-of-motion techniquei^. Analytic 
expressions for the gap energies at the F point were also 
derived assuming sufficiently small D2 , J( and J2 : hui ~ 
2S^/—D2Xi, huj2 



2Sy/-D2\i - 2{J{ 
and 



J2)Ai, huj3 



S^ 2{Di\2-D2\3-2DiJ[) 

S^2{Di\2 - D2X3 - Di{J[-A4) - 2(J{ - 4)\2) 
(from low to high energies), where Ai — 
Di + (3/2) Ji + 3J2,A2 = (3/2) Ji + 3J2 and 
A3 = 2Di + (3/2) Ji + 3J2. Fitting the calcula- 
tions to the data, we obtained Ji = —3.4(2) meV, 
J2 = -2.02(7) meV, J{ - J2 = 0.014(2) meV, 
Di = -0.28(1) meV and D2 = 0.0007(6) meV. Solid 
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FIG. 5: Constant- ?icij = 1 meV scan (a) along the (1 -|- 
h,~2h,0) direction and (b) along the (1,0,1) direction at 
T = 40 K. The horizontal bar in (a) represents the instrumen- 
tal resolution. T-dependence of (c) the intrinsic peak width, 
K, obtained from the constant-Tio; — 1 meV scans the along 
(1 -|- h, —2h, 0) direction, and (d) the relaxation rate Pqei ob- 
tained from the constant-Q ~ (1,0,0) scans shown in Fig. 3. 
Lines are explained in the text. 



lines in Fig. 4 represent the calculated dispersion 
relations with = 0- The good agreement confirms 
the validity of the model Hamiltonian. In the above, we 
could only determine the difference J( — J2 for interplane 
interactions. From the analytic expressions we see that 
huj3 and tiuji must be accurately determined in order to 
obtain J[ and separately. However, this was impos- 
sible since they appear as one peak at huj = 5.3 meV in 
Fig. 3(d) or 3(e) due to the insufficient energy resolution 
at high energies. Since the splitting between Hlo^ and 
huj4 becomes sensitive to J{ (or J2) at Q = (1.05,0,0), 
we performed a constant-Q scan at this Q and found 
an almost resolution-limited peak at huj = 5.4 meV. 
This requires the splitting to be less than the energy 
resolution AE = 0.5 nieV, and consequently an upper 
limit of 0.08 meV is obtained for J[ and Jj. Hence, the 
interplane interactions are at most 2.4 % of the inplane 
interaction Ji, evidencing the good two-dimensionality. 
One may note that Ji ^ J2, which makes YMnOs rather 
closer to the ideal TLAFM than a system of weakly 
coupled trimers. Our Ji ~ —3.4 meV is one order of 
magnitude smaller than J deduced in a recent Raman 
scattering studjii^. They obtained J ^ —140 cm^^ 
(r^ —17 meV) by assigning a broad peak appearing at 
1800 cm~^ (~ 220 meV) to two-magnon scattering. 
However, our results clearly show that the peak cannot 
be due to the two-magnon process because the band 
width of the one-magnon branch is only about 16 meV. 
Their broad peak at 220 meV must be vibrational or 
electronic in origin rather than magnetic. 

Now let us turn to the low energy quasielastic con- 
tinuum observed below Tn, clearly seen in Fig. 3(b) 
and 3(c). For T = 7 K << Tn, Fig. 3(e) shows a 
hu! = 0.55 meV mode at Q = (1,0,^) with / 7^ which 
is due to inplane transverse spin fluctuations. The in- 
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plane transverse fluctuations cannot, however, appear at 
Q = (1,0,0) because the polarization factor in Eq. (1) 
vanishes for the ordered spin structure in YMnOa. Note 
that such a mode does not show in Fig. 3(d). Therefore 
we rule out the inplane transverse spin fluctuations as the 
origin of the quasielastic continuum existing in the Neel 
phase. In order to understand the continuum, we per- 
formed constant huj = 1 meV around (1,0,0) at several 
temperatures and along different Q-directions. Shown in 
Fig. 5(a) and 5(b) are representative scans at T = 40 K. 
A nearly resolution-limited peak is seen along the inplane 
{1 + h, —2h, 0) direction whereas the intensity is inde- 
pendent of I perpendicular to the plane. These indicate 
that the quasielastic component is purely 2D in nature, 
well localized at the 2D antiferromagnetic zone center. 
Fig. 5(c) shows the temperature dependence of the in- 
trinsic peak width along the + —2h, 0) direction. For 
T > Tn the width decreases almost linearly, whereas it 
becomes nearly resolution-limited below Tn, indicating a 
large inplane correlation length at low temperatures. The 
energy width Fqei of the quasielastic peak is also shown 
in Fig. 5(d). Fqd decreases as T decreases down to Tn 
and saturates to a value of Fqei ~ 1-2(1) meV below T/v. 
It is surprising that the fast spin fluctuations with the 
characteristic time scale of Tqd = h/Tqci ~ 0.55(5) ps 
coexist with the long-lived spin-waves in the Neel phase. 

What is the origin of the fast 2D fluctuations in 
the Neel phase? Recently, a similar quasielastic peak, 
called central peak, has been found in numerical simu- 
lation studies on 2DXYTLAFMgi&. Theoretically, such 
a central peak has been commonly seen in 2D XY 
spin systems, triangular or non-triangular, and is re- 
lated to the vortex dynamics intrinsic to the KT phas^i. 
We, thus, fitted our k and Fqoi to the phenomenolog- 
ical functions: k — kq exp(— fe/^/r) and F = Fq + 



A^(A/n)e-2f/y^[(y2 - l)(ln(fcBrKT/A)/2 + b/^)]'/^, 
where r = (T - rKT)/TKT and A = JS^a^nl/A. Here 
if A = 1 and Fq = 0, k and F reduce to the analyti- 
cal expressions for the 2D XY square lattice systemic. 
The best fit (solid lines in Fig. 5(c) and 5(d)) was ob- 
tained with Tkt = 11(10) K, = 10(4), kq = 4(1) A'^, 
A = 0.07(1) and Fq = 1.2(1) meV. The fit reproduces k 
and F well for the entire temperature range, suggesting 
that the quasielastic peak is reminiscence of the vortex 
dynamics. Coexistence of the magnons and quasielastic 
peak suggests competition between the Neel phase fa- 
vored by the weak interplane interactions, and the KT 
phase intrinsic to the 2D XY spin system at low temper- 
atures. It remains to be seen whether or not the prefactor 
A being smaller than 1 and the nonzero Fq for the relax- 
ation rate are intrinsic to the 2DXYLTAFM or are due 
to the competition between the two phases 

In summary, using inelastic neutron scattering mea- 
surements on a single crystal of the hexagonal antifer- 
romagnet YMnOs, we have found in the Neel phase a 
central peak at the 2D AFM zone center which bears 
characteristics of the KT phase intrinsic to the 2D XY 
spin systems. Understanding in detail how the Neel and 
KT phases compete and change the nature of the dy- 
namic spin correlations would require further theoretical 
and experimental studies in the 2D XY spin systems. 
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